Quantum phase transitions and magnetization of an integrable spin ladder with new 

parameters in bridging to real compounds 
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We study the field-induced quantum phase transitions (QPT) and the relevant magnetic proper- 
ties of a spin- 1/2 two- leg integrable spin ladder (ISL), of which the system parameters in bridging 
to the real compounds are determined by setting the extra interactions in the Hamiltonian of the 
ISL relative to the Heisenberg spin ladder to vanish in the expectation in the ground state (GS). 
Such an ISL analytically has the correct leading terms of both the critical fields of the two QPT's 
as in the real strongly-coupled compounds: q^bH c \ — J± — Jm and <?/iB-ff C 2 = J± + 2Ju in terms 
of the experimental leg (Jii) and rung (J±) interactions. The symmetric magnetization inflection 
point is located at g^sHip = J± + Jii/2. The magnetizations for the GS and at finite tempera- 
tures, as well as the susceptibility, show good agreements in various comparisons with the finite-site 
exact diagonalization, the transfer-matrix renormalization group numerical result, the perturbation 
theory, and the compounds (5IAP) 2 CuBr4-2H 2 0, Cu2(C 5 Hi2N2)2Cl4 and (C 5 Hi2N) 2 CuBr 4 . 



I. INTRODUCTION 

Recently spin ladders [1] have attracted consider- 
able interest due to increasing experimental realiza- 
tions of ladder-structure compounds, such as SrCu 2 03 
[2], Cu 2 (C 5 H 12 N 2 ) 2 Ci4 [3], (5IAP) 2 CuBr 4 -2H 2 [4], 
(C5Hi 2 N) 2 CuBr4 [5], etc.. Among the most interest- 
ing problems related to spin ladders are the existence 
of a gap, field-induced quantum phase transitions (QPT) 
and relevant magnetic properties. As it has been ob- 
served experimentally, strongly-coupled spin ladders are 
gapped and the magnetic field H induces two QPT's 
with two critical fields. The first critical field H c \ closes 
the gap and the second one -ff c2 fully polarizes all the 
spins. In terms of the weak leg (Jy) and strong rung 
(J_i_) interactions, the leading terms of the critical fields 
axegnsHci = J_i_ — J\\ and g/j, B H c2 = J± + 2J\\ [6-8,3,9], 
which coincide with the experiments. Many spin lad- 
der compounds can be described by the standard Heisen- 
berg spin ladder (HSL), though some other compounds 
possess spin frustration and contain diagonal leg interac- 
tions. Theoretically, properties in a magnetic field have 
been discussed by numerical methods such as exact di- 
agonalization [10], transfer- matrix renormalization group 
(TMRG) technique [11], quantum Monte Carlo simula- 
tions [12], as well as other methods for strong-coupling 
limit like perturbation expansions [13], bosonization [14] 
and mapping into XXZ Heisenberg chain [9]. Recently 
some effort was exerted in bridging the integrable spin 
ladder (ISL) [15] to real spin ladder compounds by in- 
troducing a rescaling parameter 7 in the leg interaction 
[16,17], the high-temperature expansion (HTE) [18,19] 
was developed for the isotropic ISL [17] and also ex- 
tended to the anisotropic ISL with an XYZ rung interac- 
tion [20]. For a real spin ladder material, it is favorable 
to form rung states, since it is unfavorable for the weak 



leg interaction to take apart the strongly-coupled rung 
states. It happens that the ISL is based on the permu- 
tation of the formed rung states without breaking them. 
Therefore, if the system parameters are properly cho- 
sen, the ISL can provide some reasonable information, 
while from the ISL the analysis, for its characteristic fea- 
tures such as the ground state (GS), excitations and the 
gap, the QPT's, the magnetization inflection point (IP) 
and thermodynamical properties in finite temperatures, 
is quite convenient. The rescaling parameter 7 intro- 
duced in the ISL yields a gap A = J± — 4J||/7 and the 
two critical fields g^BH c \ = A, g[iFjH C 2 = J± + 4J1/7 
[16,17]. As far as the leading terms are concerned, the 
value 7 rs 4 for strong coupling fits the gap (and conse- 
quently H c i) of the HSL [16], but the second critical field 
g/is_ff c2 = Jl + Jii does not agree with the experiment- 
coinciding result giisH C 2 = J_l+2Jii [3,9]. Another point 
is that the magnetization symmetric inflection point (IP) 
Hip = (H c \ + H c2 )/2 takes the value J±/{gnB) for all 
choices of 7, the deviation from the experimental result 
Hip s» (J± + J\\/2)/(gnB) will become more consider- 
able, when the leg interaction gets stronger. A more 
appropriate way based on more physical origin for decid- 
ing the system parameters is lacking and expected to lay 
a closer bridge of the ISL to the real spin ladder com- 
pounds. 

In the present paper, we shall consider an ISL where 
the system parameters are determined by the vanish- 
ing of the expectation value of the different interac- 
tions of the ISL relative to the HSL, for the GS in 
the gapped phase and the fully-polarized phase. We 
find that such an ISL analytically has the correct lead- 
ing terms of the critical fields both at the first and 
the second QPT's: gnBHc\ = Jl — J\\ an d 3/isi? C 2 = 
Jl + 2 Jii. The magnetization IP is located at g/isHip = 
J± + Jj|/2. Our result of the magnetizations for the 
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GS and at finite temperatures compare well both with 
numerical computation from the TMRG as well as the 
finite-spin exact diagonalization and with the experi- 
mental data for the compounds (5IAP)2CuBr4-2H20, 
Cu 2 (C 5 Hi2N2)2Cl4 and (C 5 Hi2N) 2 CuBr 4 . The magnetic 
susceptibility coincides with the experimental data, and 
comparison with perturbation theory for strong coupling 
shows a good agreement. 

II. EXACT EQUIVALENCE OF HSL AND THE 

ISL FOR H > H C 2 

Many spin ladder compounds are described by the HSL 
with the Hamiltonian 

Hhsl =H\\ + H^+M, (1) 

L 

H\\ ~ ,J \\ || ' || = * ' i+1 * ' 

i=l 

W± = J± £ £ ■ A< = -g^ B H {SI + T[) , 

i i 

where S and T are spin-1/2 operators for the two legs, 
the three parts 7Y|| , 7i± and M. denote the leg interaction, 
rung interaction and the Zeeman energy, respectively. Jy 
and J± are the leg and rung coupling strengths. There 
are totally L rungs and H is the external magnetic field. 

Before considering for all values of the magnetic field, 
we shall first take a simple look at the fully-polarized 
phase and show that the HSL is equivalent to an ISL 
for H > H C 2- For convenience we label the rung sin- 
glet as |0i) = |t| - IT), /V2 and the rung triplet as 
\U) = 1% \2i) = III + 11), /V2, 13;) = HI), where 
the subscript i denotes the i'th rung. The rung interac- 
tion leads to an energy — 3Ji/4 for the singlet and Jj_/4 
for the triplet. After H c2 all spin are polarized so the 
ground state is ^fp = | I1I2 • ■ ■ 1l) for a ladder with to- 
tal L rungs, which is exactly the eigenstate of the HSL 
Hhsl- 

When the field is lowered down to the critical point H c2 
from the fully-polarized phase, spins begin to flip down. 
More elementally the spin flips are linear combinations of 
the strongly-coupled rung states. Less excitations lead to 
lower energy, so one the-most-favorable state will enter 
the GS first, which indicates the occurrence of the QPT. 
Since the rung state |3) has higher energy in the pres- 
ence of the field, we only need to consider |0) and |2) for 
the critical point. The lowest excitation can be trapped 
out from a subspace where there are at most one |0) and 
|2), all possible neighbors collected are {| lili + i), |ljOj + i), 
|L2 l+1 ), \0A +1 ) (and |0il <+ i), |2<l i+1 >, |2jO i+ i»}. For 
all these neighbors we exactly have an equivalent relation 

ft|l' t+1 = -Pi.i+i, ( 2 ) 

where P M+i = [2S l ■ S i+1 + \){2Ti ■ f i+1 + |) is the per- 
mutation operator. P»,j+i exchanges the rung states <pf. 



Pis+i I fif'i+i) —\ fifi+i)- Consequently the Hamilto- 
nian in this subspace can be rewritten as 

Hhsl - ]T f^^.i+i + J±& ■ % - 9VbH(S* + T*)\ 

i ^ ' 

which is actually an ISL [15]. As there is at most one |0) 
and one |2), the exact eigenenergy from the Bethe ansatz 
approach [24] is 

E = E FP + {gn B H) (No + N 2 ) 

N +N 2 

-o J ll E 72-T- J ^o, (3) 
1 j=l Pj + Z 

where E FP = L (|Jm + jJ± — g^sH) is the fully- 
polarized eigenenergy, No and N 2 (taking the value 
or 1) are the numbers of states |0) and |2), respectively. 

Apparently both |0) and |2) are gapful if gfisH > 
max{Jj_ + 2J||, 2Jy}, and the GS is fully-polarized. The 
value /j, = corresponds to the critical field <7/xsi? C 2 = 
max{Jj_ + 2J||,2J|i}, which reaches the same result de- 
duced from the instability of spin wave spectrum in the 
ferromagnetic phase [3] and that obtained from mapping 
onto an XXZ Heisenberg chain [9] . The expression of the 
saturation point is exact [3], more exactly speaking, for 
an even number of rungs. The one-particle Bethe ansatz 
wavefunction corresponding to fi = provides 

L 

*c2 = ^(-ir|lll2---6 i ---U), (4) 

i=l 

where 6, = 0, 2 according to favorable energy. For an- 
tiferromagnetic spin ladder with positive J± and Jm, |2) 
has higher excitation energy, then hi is the singlet |0) and 
g/j.E-H C 2 = J±_ + 2Jii. The \I/ C 2 will be used in our fur- 
ther discussion. For an odd number of periodic rungs, 
the term 2Jm in the critical expression is replaced by 
(1 — cos[(l — 1/L)7r]) Jm, as /j, = is not the solution in 
periodic boundary conditions. But this size effect will 
completely disappear in thermodynamical limit. 



III. ISL IN BRIDGING TO THE COMPOUNDS 
IN ALL FIELD VALUES 

After the simple look at the fully-polarized phase in the 
previous section, we proceed to discuss all values of the 
field applied on an antiferromagnetic spin ladder. The re- 
lation 7Yj|' 1+1 = \Pi t i+\ from (2) docs not hold for neigh- 
bors |0j0j+i), |2j2j + i) or |lj3j+i), though it is exact for 
all the other kinds of neighbor states. To find a model 
effective for all value of fields, we shall consider such an 
ISL with the difference in Ho 

H ISL =Ho+H±+M, (5) 

nj j \ "> nji,i+l n/i,i + l . nji,i+l 

Ho - J\\ 2^ n o > n o - + H extra , 
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where 



nji,i+l 

^■extra 



^i±l 4+1 + |$ • ?i + SWi • f l+1 ) + c. 

(6) 



The rung interaction part TL±_ and Zeeman term M. arc 



exactly the same as the HSL Hhsl- If ti extra — > t^ 10 
H'q 1+1 will reduce to the leg interaction of the HSL. The 



extra part W^ct™ makes the model integrable, but the 
biquadratic exchange Si ■ Si + \Ti ■ T i+ i [21] in P,,i+i docs 
not appear in the HSL, though it is as weak as the leg 
interaction. We set and balance the parameters 7, a and 
c to reduce its influence. To minimize its extra effect, we 
consider in such a way that it vanishes in expectation 



extra/ 



0. 



(7) 



i=i 



as much as possible in the GS. Our consideration aims 
to make the model exactly soluble on one hand, while on 
the other hand, the extra effect is minimized as much as 
possible. 

More complicated balance terms such as S?Tf, S^ + T[ 
can be taken into account in but they break the 

isotropy. Here we consider the simplest case (6) which 
preserves the isotropy when it is giving quite good re- 
sults. To see the integrability, the Hamiltonian can be 
rewritten as 

Hisl = J\\ ^C^ 1 + c) + (Jl + aJ\\)^2 St-Tt+M 

i ^ i 

(8) 

where Pi,%+\ forms the permutation operator in the SU(4) 
integrable model [22,23] and M. is the Zeeman term as in 
the HSL. The model (8) essentially is an ISL [15] which 
can be exactly solved by Bethe ansatz approach [24] . The 
exact eigenenergy reads 



M 



(i) 



' j Mj + 4 fe=o 

3 1 
E = --(J± + aJy), Ei = -(Ji + aJ\\) - g/j, B H, 

E 2 = ^(Jj_ + aJ||), £ 3 = ^{J± +aJ\\)+giJLBH, 

where Ek is the single-rung energy of the singlet (Eo) 
and triplet (E\, E2, E3), Nk is the total rung number 
in the corresponding rung states. The rung state with 
lowest Ek is chosen as the reference state, is the 

total number of rungs occupied by the other rung states. 
The singlet is chosen as the reference state before the IP 
g^BHjp = (J_l + atJ\\) (decided by E = E\), while |1) 
is the reference state thereafter. The gapped phase lies 
in the former case and the fully-polarized phase locates 
in the latter. 



The gapped phase of the ISL is composed of the singlet. 
The magnetic field lowers the energy of the triplet com- 
ponent |1) which enters the GS when the gap is closed at 
the first QPT. Further increase of the field brings all the 
components of the singlet out of the GS and another gap 
opens at the second QPT (we refer the gapped phase to 
be at H < H c i, though the fully-polarized phase is also 
gapped). Since both the gapped phase and the fully- 
polarized phase only consist of one rung state (|0) or |1), 
respectively), which is chosen as the reference state, 
in these two phases corresponds to the number of the ex- 
citations to the other three components. Then from the 
eigenenergy (9), it is very easy to find the excitation gap 
and the two critical fields 



gfJsHci = (Ji + aJ||) - -J|| 

4 

gHBH c2 = (Ji + aJ||) + - J|| 



(10) 



which, respectively, closes the gap of the gapped phase 
and open the gap of the fully-polarized phase. The 4/7 
terms come from the lowest energy in the excitation 



band, \i 



(i) 



corresponds to the bottom of the en- 



ergy band. The critical fields as well as the IP realize the 
relation Hjp = i (H c i + H C 2), which is the magnetiza- 
tion symmetric point and agrees with the experiments, 
as discussed in the following. 

Experimentally strongly- coupled spin ladder com- 
pounds have two similar QPT's, separating the applica- 
tion of the external magnetic field into three phases. As 
usually described by the HSL, the fully-polarized state 
>J>p exactly is the system GS after the second phase 
transition, while the singlet is the most favorable state 
and dominates overwhelmingly in the gapped phase be- 
fore the occurrence of the first QPT for the antiferromag- 
netically strongly-coupled real spin ladder compounds. 

Among the total three phases, the extra-interaction ef- 
fect of the ISL can be minimized explicitly in two phases, 
gapped phase and the fully-polarized phase. Fulfilling the 
relation (7) in these two phase gives rise to 



- - I + -a + c = 0, H>H c2 ; 
7 2 4 



1 3 

7 4 



a + c = 0, H < H cl , 



(11) 
(12) 



where the term 1/2 in the former equation comes from 
the expectation of the HSL leg interaction, while in the 
latter equation this expectation vanishes in the singlet 
phase. The solution is direct to obtain and we yield to 
the limitations 



1 3 1 

78 2 



(13) 



Under these parameters, the GS has the exact equiva- 
lence Hhsl = Hisl for the fully-polarized phase H > 
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H C 2 and the expectation equivalence (Hhsl) — (ft-isi,) 
for the gapped singlet phase H < H c \. Then the IP is 
determined from the definite a in (13) 

gH B H IP = Jj_ + , (14) 

which, as it turns out in the following comparisons, pro- 
duces well the experimental value and TMRG numerical 
result. Note here that this experiment-coinciding IP ex- 
pression is an analytic result from our consideration, and 
instead of from empirical parameter choosing, the pres- 
ence of the term i Jy comes from the only requirement on 
the extra-interaction-effect minimizing (7) of the above 
two phases. We shall decide more definitely the parame- 
ter 7 by considering the critical points of the QPT's. 

A. Equivalence at the 2nd QPT 

Since the exact analysis for the HSL is available at 
the second QPT as discussed in the previous section, we 
first consider such a case that the extra interactions of 
the ISL also vanish at the second critical point, in addi- 
tion to the requirement on the gapped phase and fully- 
polarized phase. Although ft'gxtra does not vanish for 
a single neighbor-site pairs |ljOi + i) at H c2 , the action 
of these extra terms together can cancel each other. In 
fact, noticing the relation (2) for the rung-state neigh- 
bors |lil i+ i) and | liOi + i) involved at the second critical 
point, one can easily obtain 

X)K££«*c2 = [(\-\ + J )(L 4) + c£]* c2 . (15) 
i=i ~ 

Setting the extra-interaction effect to vanish 
J2i=i % extra^ c2 = 0, in addition to the limitations (13) 
from the afore-discussed two phases, leads us to the 
definite parameters 

8 1 

7= -, <*=-, c = 0. (16) 

Thus the relation (2) is replaced by TC^ +1 = H.Q l+1 for 

fully-polarized phase H > H C 2 and J2 i HV l+ = 

J2ift(> 1+1 at the critical point H c2 . Such a simple con- 
sideration immediately gives the two critical fields 

g/i B H c i = Jj_ — Jy, 

gHBH c2 = J± + 2J ll . (17) 

As it is expected, H C 2 is the same exact result as we dis- 
cuss in the previous section. Note that here H c \ for the 
first QPT is actually composed of the leading terms from 
the perturbation theory and has been often used in ex- 
periments, while we have not imposed any limitation on 
the first QPT and it is a by-product of the consideration 
of the second QPT. 



The magnetization can be obtained by thcrmodynam- 
ical Bethe ansatz (TBA) [25,26] which is also very conve- 
nient for analysis of the field-induced QPT's [16,27,20]. 
For the ISL under consideration, the magnetization be- 
tween the two critical fields in zero-temperature is de- 
cided by the following TBA equations of the dressed en- 
ergy and density, 

e« = 5 «-2^ ai -a 2 * e ( 1 )-, 

=ai-a 2 *P (1) ", (18) 

where gW = ±^[(l+a)J±-gn B H\/J\\ for ±(H IP -H) > 
0, a n (fi) = [n/([i 2 + n 2 /4)]/(27r), and the symbol * 
denotes the convolution. The Fermi sea with negative 
dressed energy e' 1 ^ forms the system GS. Here the GS 
only involves one branch of dressed energy, since the GS 
of the ISL is composed of the singlet in the absence of the 
field, while only the component |1) is lowered in energy 
level and brought down to the GS by the applied field. 
Correspondingly the gap of the excitation to the compo- 
nent 1 2) and |3) never closes during the gapless competi- 
tion between |0) and |1) in the GS. We present an exam- 
ple of the GS magnetization in Fig.l. The symbols are 
data figured out from Ref. [10] in the study by Lanczos 
exact diagonalization on 2x12 and 2x16 HSL with peri- 
odic boundary conditions, the ratio of the leg and rung 
interaction strength is J\\/J±_ = 0.2. The solid curve is 
our result which coincides with the finite-site exact di- 
agonalization. There is some small discrepancy at the 
first QPT due to existence of higher-order terms in H c i 
(the discrepancy will be reduced when the higher-order 
terms are picked up, as one will see from GS magnetiza- 
tion compared with T=0.02 TMRG result in Fig.5 and 
Fig. 7). As our consideration provides the same H c \ and 
H c2 as the HSL in leading terms, the coincidence of our 
result with the experiments is also expected. Two exam- 
ples are given in Fig.2A and B, respectively for the com- 
pound Cu2(C 5 Hi2N 2 ) 2 Cl4 [3] with ratio J\\/J±_ w 0.18 
and (C 5 Hi2N) 2 CuBr4 [5] with J {l /J ± w 0.29. 

The high-temperature expansion (HTE) [18,19] has 
been developed for the isotropic ISL [17]. Then the 
temperature-dependent free energy per spin takes the 
form 

OO j 

n=0 ' 

where the expansion is carried out for per rung and the 
factor 1/2 is added for per spin, k B is the Boltzmann 
constant. We present some orders of coefficients, 

2Q 3Q 6Q 2 3Q_ 

C -ln Q+ , Cl -Ql> C2 -Ql~Ql + ${' 

10Q_18Q^ 80Q^ 8Q_ 24QQ_ 
° 3 ~ 3Q 2 + Q\ + 3QI + Q\ Q\ + Q| ' 
35Q 205Q 2 120Q 3 140Q 4 55Q_ 
4 " 12Q^ " ~6Q\ + ~Ql Q\~ + ~4Ql 
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100QQ- 180Q 2 Q_ AhQ 2 __ 15 40Q 



Q+ 



Q+ 



2Q% Q% Q% ' 



A = J±- J\\+dJ h d = a(^-) + b(^-) 2 , 



(20) 



where 



Q = 2cosh(i/?J) +4cosh(i/3J)cosh(/3/i), 
g± = 2e( ±/3j / 4 ) cosh(i/? t 7) + 2 e fr/"/ 4 ) C osh(/3/i), 

/i = g/J, B H, J = (J_l + «J||) and /3 = l/(k B T). 
One can get higher orders for lower temperatures. 
We compare our result of (16) with the compound 
(5IAP) 2 CuBr 4 -2H 2 [4], with ratio J^/J ± w 0.077 ± 
0.015, the magnetizations M = —dF/dH and the mag- 
netic susceptibility \ = —d 2 F/dH 2 are presented in 
Fig. 3 and Fig. 4, respectively. The comparison shows a 
good agreement with the experiments. The unit of the 
susceptibility in comparing the experiment is determined 
by the Curie constant C, which can be clearly seen when 
we expand the susceptibility in (27). The Curie con- 
stant is decided, usually for molar susceptibility, when 
the g factor is measured by the experiment. For the com- 
pound Cu2(C5Hi2N2)2Ci4 [3], we have shown in Fig.2B 
the coincidence of our GS magnetization with the exper- 
imental magnetization of the lowest available tempera- 
ture. This compound possesses nonlinear torque mag- 
netization [29] and possible weak diagonal interaction 
[10], which may be the possible reason that causes the 
observed asymmetric magnetizations at different higher 
temperatures. So here we shall not compare the magne- 
tizations for this compound at higher temperatures. The 
compound (C5Hi 2 N) 2 CuBr4 [5] has stronger leg interac- 
tion ( J\\/J± ~ 0.3 ), the higher-order terms of the gap 
become more considerable. The small difference resulted 
from the high order terms can be seen from Fig.2B, the 
dotted line and the solid one are the magnetizations in 
the absence and presence of the high-order terms, respec- 
tively. In the following we shall pick up the high-order 
terms. 



B. Fitting of the 1st QPT 

Up to now our consideration has exhibited the correct 
leading terms of the critical fields. When the leg interac- 
tion is stronger, the higher-order terms in the gap (equal 
to g^sHd) become more considerable. On the other 
hand, the fluctuation (H 2 HSL ) — {Hhsl) 2 in the singlet 
phase is very large, treatment of the first QPT may in- 
fluence the system properties more sensitively than the 
second QPT. Unlike in the second QPT, the exact solu- 
tion of the HSL is not available in the first QPT, requiring 
Ti-extra^ci = is not explicitly executable. So we di- 
rectly adjust the parameters to fit the gap, which is also 
a result of (Hhsl) — (Hisl)- The higher orders terms in 
the gap of the HSL converge slowly in the perturbation 
theory, but it was found [12] that a two-parameter gap 
A fits the numerical data of Greven et al. [30] very well 



where a = 0.6878, b = -0.1861. The dJ\\ comes from 
the high-order terms in the gap and here we have used a 
label d for our convenience. Setting 



1 _ 3 d _ 1 

7 ~ 8 ~ 4' a ~ 2' 



d 

C= 4' 



(21) 



will give rise to this gap at the same time when the re- 
quirement in (13) of the ISL is fulfilled. In Fig. 5 we 
compare the magnetizations with the TMRG numerical 
result [11] in various temperatures, for Jn = 1, Jj_ = 5.28. 
The solid lines are plotted according to the correspond- 
ing temperatures of the TMRG data, the broken line is 
plotted at T= compared with the lowest temperature 
T=0.02 TMRG result. As one can see, the agreement is 
quite good for all the temperatures. For the comparison 
with the experiment, Fig. 6 presents the magnetizations 
for the spin ladder compound (CsH^N^CuB^ [5] with a 
stronger leg interaction (Jn/J± w 0.29, d w 0.18). As the 
HTE is not valid for 0.7K, we also calculate the magne- 
tization of GS (0.0K) by TBA, a more clear comparison 
has been presented in 2B. The figure demonstrates a good 
agreement of our result with the experiment. The sus- 
ceptibility in a low field is given in the inset, our result in 
the solid line coincides with the experimental observation 
except for some small deviation at the hump. 

The interesting IP is an invariant point of the 
magnetizations at low temperatures. The compound 
(5IAP) 2 CuBr 4 -2H 2 [4] in Fig.3 indicates the IP, though 
it suffers from the impurities. As the compound 
(C5Hi2N)2CuBr4 is free of the impurity or diagonal in- 
teraction, it exhibits an unambiguous IP at the half- 
saturation. The TMRG numerical calculation for the 
theoretically ideal HSL also shows the IP. From the ISL, 
it is direct to understand the IP [20,16,17]. Two reasons 
are both necessary for the IP: (i) equal energy of the only 
two components competing in the GS; (ii) large excita- 
tion gap relative to the temperatures. Actually only two 
components, the singlet |0) and the triplet component |1), 
compete in the GS, as we mentioned below the dressed 
energy equation. The IP locates at the point where the 
reference state in an increasing field converts from |0) to 
|1). At this point the two components have the same 
rung energy Eq = E\ and same proportion No = N\ in 
the GS. The excitations to the other components |2) to 
1 3) are gapful, the gap are exactly available [20] 



A /P = (J J _+aJ||)-(2hi2)J||/7, 



(22) 



which is large due to the strong rung coupling J±. Low 
temperatures can hardly stimulate excitations to |2) or 
|3), while each of |0) and |1) with the same energy still 
occupies half the total rungs. As a result, the IP is invari- 
ant with a half-saturation magnetization under low tem- 
peratures. Higher temperatures will excite |2) as well as 
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1 3) and lower the magnetization, consequently the mag- 
netization curves do not go through the IP, as one can 
see partly from the temperature T=2.5 and more clearly 
from T= 5.0, 10 in Fig.5. The location of the IP is well 
produced by our analytic expression (14) for the com- 
pound and the TMRG data. For the unambiguous IP of 
the compound (C 5 Hi2N) 2 CuBr4 (g = 2.13, J± = 13.3K, 
J|| = 3.8K) our result Hjp = 10. 6T is identical with the 
experimental observation; the result for J± =5.28, J|| =1 
from TMRG keeping 200 optimal states provides an IP 
at H IP w 5.83 (H cl w 4.38, H c2 w 7.28 and g/j, B is in- 
corporated in the field as in Rcf. [11]) when our result is 
Hjp = 5.78 with a very minor discrepancy about 0.9%. 



some small improvement at the peak of the low-field sus- 
ceptibility (given in the inset of the Fig. 6). But for the 
theoretically ideal HSL under the couplings J± = 5.28, 
J|| = 1, the value d w 0.124 from (20) gives the IP 
g^pHjp = 5.84, which minimizes the discrepancy from 
the TMRG result gfipHip = 5.83 [11]. As our previous 
considerations with simple parameters have given quite 
good results, here we just give a comparison for the GS 
magnetization with the TMRG data in Fig. 7. 



D. Susceptibility compared with perturbation 
theory of HSL 



C. Fitting both the two critical fields 

Under the parameters (21) for fitting the gap, the 
extra interactions at the second critical field will be 
J\\J2i=i^extra 1 ^'c2 = — dJ|i* C 2j which does not vanish 
though it is small. Nevertheless, this can be also resolved 
by more consideration. Previously we only take the dom- 
inant singlet component into account in the gapped phase 
of the HSL, as a simple physical picture as in the ISL. 
When the HSL-type leg interaction becomes less weak, 
small amount of other components will be mixed into the 
gapped singlet phase of the strong coupling limit. There- 
fore the average per-leg expectation value of the HSL leg 
interaction does not vanish as in a pure singlet phase. 
This contributes to another modification parameter S. 
As we are using the ISL to describe the compounds, the 
(12) becomes 
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-a + c = 5, H < H c i, 



(23) 



instead. In limitation (12), S is null in our previous 
simpler consideration that singlet component dominates 
overwhelmingly in the strongly-coupled HSL. Combining 
(23) in the gapped phase, the gap fitting g^pHd = A 
and the extra- interaction minimizing (11) in the fully- 
polarized phase, one gets 



1 

7 



d+5 1 S d . , 

— «=2"*' c =2 + 4- (24) 



If setting S = —d/2, one will regain the vanishing of 

the extra-interaction J2f =1 %extra^c2 = for the second 
critical field. In such case we will have the critical fields 



g^ B H c i = Jj_ - J|| + dJ\\, gfiBH c2 = Jj_ + 2Jy, (25) 
as well as a new IP 



I d 
gn B Hip = J± + -J|| + 



(26) 



For the compound (CsHi 2 N) 2 CuBr4 [5], these results 
does not give improvement to the magnetizations except 



In bridging to the experimental data, the unit is a prob- 
lem needing some special attention. The unit of the sus- 
ceptibility can be clarified if the expansion is carried out. 
One can obtain the expansions of the susceptibility for N 
spins easily from the free energy expression, the lowest 
two orders take the form 



C 
T 



(27) 



where C = N g 2 n 2 B / (Akp) is the Curie constant. In the 
expression (27) the part in the square brackets is dimcn- 
sionless. The Curie constant is provided by the experi- 
ments corresponding to a measured value of g, e.g., C = 
0.412 emu K/mol for the compound (5IAP) 2 CuBr 4 -2H 2 
(g = 2.1) [4]. As one knows, the lowest order \ — C/T is 
the Curie law, which guarantees the coincidence with the 
experimental data in high temperatures (T — > 100K for 
the compounds discussed in the present paper) . Compar- 
ing (27) with the Curie- Weiss Law \ = C/(T—Q) implies 
a Curie- Weiss temperature Qisl = — [J± + (2/j+a) J||]/4 
for the ISL (8). The HSL has a Curie- Weiss temper- 
ature Qhsl = —(J± + 2J||)/4 [12]. Our parameters 
7 = 8/3, a = 1/2 provide ISL = -(J± + 5J||/4)/4, 
which is close to that of the HSL. The difference of Qisl 
from &hsl will be more ignorable when the leg inter- 
action Jii is relatively small in a strongly-coupled spin 
ladder, which coincides with the good agreement of sus- 
ceptibility for the compound (5IAP) 2 CuBr4-2H 2 0. 

For the strongly- coupled HSL, it was found that the 
expression of zero-field susceptibility derived from per- 
turbation theory [12,13], 



X 



^4Cf 1 Jj 
T + e B J ± 



IB 



(3 + e B ) 2 



h 
J± 

j± 



2 r 



3 r 



4(3 + e B ) 3 
3B(e 2B - 1) 



8(3 + e B f 
9B 2 e B (l + 3e B )-B 3 (7e 2B 



9e B - 12) 



12(3 + e B ) 4 



(28) 
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reproduces with high accuracy the quantum Monte Carlo 
simulations [12], as long as the ratio J\\/J± does not ex- 
ceed 0.1. The definition B = Jj_/(fcsT) is utilized in 
(28). We compare our result of the susceptibility with 
(28) of the HSL in Fig. 8. For more clarity of unit, we 
use the susceptibility scaled by the Curie constant. The 
comparison for J\\/J± = 0.1 shows that our result is in a 
very good agreement with that of HSL, the coincidence 
will be even better for J\\/J± < 0.1. 



IV. SUMMARY 

We have studied the quantum phase transitions, the 
magnetizations and susceptibility of an intcgrable spin 
ladder (ISL) by proposing a new way in deciding the 
system parameters to build closer connection of the ISL 
with the real compounds. The system parameters are 
chosen to make the extra interactions of the ISL relative 
to the Heisenberg spin ladder (HSL) vanish in the ex- 
pectation of the ground state, especially for the gapped 
phase and fully-polarized phase. The analytic critical 
fields of such an ISL have the same leading terms as 
those of the experimental observations and the HSL: 
gfiBHd w J± — J|| and g\iBHc2 ~ J± + 2J||, accom- 
panied by the symmetric magnetization inflection point 
with the location analytically at q^ibUip — J± + J\\/2 
(another result g^sHip = J± + J\\/2 + dJ\\/2). A vari- 
ety of comparisons have been made with other numerical 
results and the experiments. The magnetizations coin- 
cide with the finite-site exact diagonalization of the HSL 
for the GS, and with the TMRG numerical result for 
full temperatures. For experiments, the magnetizations 
agree well with the compounds (5IAP)2CuBr4-2H20, 
Cu 2 (C5Hi 2 N 2 )2Cl4 and (C 5 Hi2N) 2 CuBr4, with the leg 
and rung coupling ratio J\\/J± around 0.1, 0.2 and 0.3, re- 
spectively. The susceptibility with the clarified unit also 
coincides with the strongly-coupled compound as well as 
perturbation theory. The coincidence in the magnetiza- 
tion and susceptibility also makes it possible to discuss 
the specific heat for the whole process of the field appli- 
cation, we shall present the discussion in another paper. 
Our method may also be helpful for other integrable mod- 
els in bridging to the experiments and real compounds. 
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Scaled field H/A 

FIG. 1. Magnetization versus the scaled magnetic field for the ground state (GS) (T=0.0K), J\\/J± = 0.2. The symbols are 
result of Lanczos exact diagonalization on 2x12 and 2x16 Heisenberg spin ladder (HSL) with periodic boundary conditions, 
the data are figured out from Ref. [10]. A is the gap of the HSL [28]. Our result (solid curve) from 7 = 8/3 and a — 1/2 
coincides with that of the exact diagonalization, there is some small discrepancy at the first QPT due to existence of higher 
orders terms in H&. The dashed lines present result of some choices from an adjustable 7 and a = 0, they all go through a 
same half-saturation point away from the exact diagonalization data. 
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Magnetic field H [T] H [T] 

FIG. 2. (A). Groundstate magnetization versus the magnetic field in comparison with the compound Cu2(C5Hi2N2)2Cl4 
[3], g z = 2.03, J± = 13. 2K, Ju = 2.4K, Ju/J± = 0.18. The symbols are experimental magnetization for the lowest available 
temperature (0.42K). The solid line is the groundstate (T=0.0K) magnetization calculated from our consideration with 7 = 8/3 
and a = 1/2. (B). Groundstate magnetization compared with another compound (C5Hi2N)2CuBr4 [5], g = 2.13, J± = 13. 3K, 
J|l = 3.8K, J\\/J± = 0.29. The temperature of the experimental data is 0.7K. The dotted line is the result of 7 = 8/3 and 
a = 1/2, while the solid line is the result including the high order terms in the gap, I/7 = 3/8 — d/4 and a = 1/2. The 
parameter d is given by (20) under certain values of J± and Jm . 
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Magnetic field H [T] 

FIG. 3. Magnetizations versus the magnetic field at different temperatures for the compounds (5IAP)2CuBr4-2H20, g = 2.1, 
J± — 13. OK, J|| = 1.0K. The data in symbols are the experimental result [4], the lines are plotted from our theoretical calculation 
(7 = 8/3 and a = 1/2). The compound contains a considerable amount of impurities, which can be seen obviously in the 
magnetization of T= 0.4K case. The impurity effect is also reflected partly in the T=1.59K case but overwhelmed by higher 
temperature in T=4.35K. There is some uncertainty of the interactions (Ju/J± = (7.7 ± 1.5)%) due to the impurity. As the 
leg interaction Jy is rather small, the relative uncertainty is quite large. Here we use the values provided in the summary of 
Ref.[4]. 
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FIG. 4. Magnetic susceptibility versus the temperature for the compounds (5IAP)2CuBr4-2H20, g = 2.1, the Curie constant 
C — 0.412 emu K/mol, J± = 13. OK, Ju = 1.0K. The symbols are the experimental observation, the solid lines are from our 
calculation. The results for 7 = 8/3 and I/7 = 3/8 — d/4 do not differ, since the higher orders terms in the gap can be 
completely ignored due to the weak leg coupling (Ju/J± ~ 0.077). Impurity contribution was subtracted in the inset, so the 
experimental curve is lowered slightly. 








FIG. 5. Magnetizations versus the magnetic field at different temperatures in comparison with the transfer-matrix renormal- 
ization group (TMRG) numerical result, J± =5.28, Jy =1 and A = 4.382. In this figure the g factor and Bohr magneton \ib 
are incorporated in the field H . The data in symbols are TMRG numerical result [11]. The solid lines and the broken line are 
plotted from our result for 1/7 = 3/8 — d/4 and a = 1/2. The GS magnetization is plotted in the broken line in comparison 
with T=0.02 TMRG data. 
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FIG. 6. Magnetizations versus the magnetic field at different temperatures for the compound (CsH^N^CuBr,!, g = 2.13, 
J± = 13. 3K, = 3.8K. The symbols denote the experimental data [5]. The solid lines are plotted from our result for 
I/7 = 3/8 — d/4 and a = 1/2. The difference for the GS magnetization arising from the higher-order terms in the gap can be 
found in Fig.2B. The inset shows the susceptibility against the temperatures in a low field 0.1T, the Curie constant C = 0.425 
emu K/mol (g = 2.13) [5] .The solid line (I/7 = 3/8 - (d + 5)/4, a = 1/2 - 5, 5 = -d/2) and dotted line (I/7 = 3/8 - d/4 and 
a = 1/2) have only small difference at the peak. 




Magnetic field H/A 

FIG. 7. Groundstate magnetizations in different parameters in comparison with the T=0.02 TMRG result [11] Jj_ =5.28, 
J11 =1 and A = 4.382. Our reults are plotted according to (13), (21) and (24), respectively, and S = —d/2 for the solid line. To 
show their small difference more clearly, we have focused the field on the region between 0.8 and 2.0. 



10 



0.08 




10 20 30 40 



Temperature T[K] 

FIG. 8. Scaled magnetic susceptibility against the temperature in comparison with the Heisenberg spin ladder (HSL), 
J\\/J± = 0.1, J± — 13K. Since all compounds used in the figures have rung interaction around 13K, we use this value here. 
For more clarity of the unit, we scale the susceptibility by the Curie constant C. The HSL curve is plotted from the expression 
(28), which reproduces accurately the quantum Monte Carlo simulations [12] for J\\/J± < 0.1. Here the high-order terms in the 
gap make little difference in our result due to the small ratio Jy/Jx. The agreement for our result with the HSL for stronger 
coupling J\\/J± < 0.1 is even better. The broken line is the result from the choice 7 = 4, a = 0, for comparison. 
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